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ABSTRACT 

Gamma rays from extragalactic sources are attenuated by pair-production interactions 
with diffuse photons of the extragalactic background light (EBL). Gamma-ray bursts 
(GRBs) are a source of high-redshift photons above 10 GeV, and could be therefore 
useful as a probe of the evolving UV background radiation. In this paper, we develop 
a simple phenomenological model for the number and redshift distribution of gamma- 
ray bursts that can be seen at GeV energies with the Fermi satellite and MAGIC 
atmospheric Cherenkov telescope. We estimate the observed number of gamma rays 
per year, and show how this result is modified by considering interactions with different 
realizations of the evolving EBL. We also discuss the bright Fermi GRB 080916C in 
the context of this model. We find that the LAT on Fermi can be expected to see a 
small number of photons above 10 GeV each year from distant GRBs. Annual results 
for ground-based instruments like MAGIC are highly variable due to the low duty 
cycle and sky coverage of the telescope. However, successfully viewing a bright or 
intermediate GRB from the ground could provide hundreds of photons from high 
redshift, which would almost certainly be extremely useful in constraining both GRB 
physics and the high-redshift EBL. 
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1 INTRODUCTION 

High-energy gamma rays traveling cosmological distances 
can interact with the diffuse UV, optical, and IR radi- 
ation fields in electron-positron pair production interac- 
tions, leading to an effective optical dep th for VHE sources 
|Nikishovl 19621 : iGould fc SchredeJl967l) . A recent technique 
in gamma-ray astronomy has been to use observations of 
relatively nearby blazars observed by Imaging Atmospheric 
Cherenkov Telescopes ( IACTs) to constrain t he EBL (e.g. 
Aharonian et all 120061 ; iMazin fc Raud 120071 : lAlbert et all 
2008ah . 

While blazars have been the primary target of efforts to 
detect the effects of EBL attenuation in high-energy spec- 
tra, another exciting possibility is to see these same effects in 
observations of GRBs. Until recently, limitations on the ef- 
fective areas and energy ranges of high-energy experiments, 
and their ability to respond sufficiently quickly to tran- 
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sient events have hindered observations at GeV energies that 
might reveal evidence of this phenomenon. 

The EGRET experiment on the Compton Gamma-ray 
Observatory (CGRO) operated with energy range 20 MeV- 
30 GeV and effective area ~1000 cm 2 . This experiment de- 
tected a total of 5 bursts above 30 MeV in 4 years of oper- 
ation , including 4 individual photons above 1 GeV l|Dingusl 
ll995T ). While these EGRET detections do suggest the pres- 
ence of a very-high energy component in the spectrum of 
some GRBs, it is difficult to draw more conclusions due 
to the small number of high-energy gamma rays detected. 
At the same time, the BATSE instrument on CGRO de- 
tected t housands of GRBs at energies between 20 keV and 
2 MeV l|Paciesas et al.lll999T ). The Swift mission has been 
finding bursts at a rate of about 8 per month since its 
launc h in December 2004 at energies between 15 and 150 
keV (jSakamoto et al.ll2008l ). 

Until the launch of the Fermi Gamma-ray Space Tele- 
scope in 2008, the only other possibility to view high-energy 
emission from GRBs was with ground-based experiments 
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- IACTs such as H.E.S.S., MAGIC, or VERITAS, or air 
shower arrays such as Milagro. While IACTs have the ad- 
vantage of much larger effective collecting areas than satel- 
lite detectors, they are limited by their low duty cycles and 
small fields of view, which make a serendipitous GRB detec- 
tion very unlikely. Therefore, these IACTs must be alerted 
to a burst event by another detector, usually a gamma- or 
x-ray satellite, and slew to its position. This introduces a 
new technical limitation on ground-based observations; the 
delay time in receiving an alert and moving the telescope 
means that much or all of the primary emission of a burst 
can be missed. Air shower arrays like Milagro do not have 
slewing issues, but are generally not as sensitive as their 
IACT counterparts, particularly at lower energies. 

Despite the difficulties involved, followup observations 
have been made of many GRBs by all major experiments 
mentioned above. MAGIC responded to 35 burst alerts 
between January 2005 and June 2008, a rate of about 
1 per month, with an average slew time of 45 seconds 
i|Garczarczvk et alJliooj . see also I Albert et al.ll2007h . These 
attempts were only able to place upper limits on the flux. 
Negative results were also found in 32 o bservations over 4 
years by H.E.S.S. l|Aharonian et alj|2009h . Most of these ob- 
servations did not begin until several hours after the initial 
detection of the event. V ERITAS has reported limits for 
a small number of bursts l|Horanl l2008). Air shower arrays 
are generally less sensitive than IACTs, but have the ad- 
vantage of much larger fields of view and duty cycles. The 
Milagro prototype, Milagrito, claimed a possible detection 
of pro mpt emission from GRB 970417A at > 650 GeV en- 
ergies (lAtkins et alj|2003j ). The redshift of this burst is not 
known, but to be detected at these energies it would have 
had to have been quite nearby, as the universe becomes op- 
tically thick due to EBL attenuation at low redshift (z < 0.2 
for our lowest EBL model) for E-y > 650 GeV. 

In this work, we focus on two telescopes that are well- 
suited to detect GeV emission from GRBs: Fermi and 
the Major Atmospheric Gamma-ray Imaging Cherenkov 
Telescope (MAGIC). Fermi contains two instruments. The 
GLAST Burst Monitor (GBM) is designed for finding 
prompt emission from GRB from 10 keV to 30 MeV, and 
the Large Area Telescope (LAT) views gamma rays from 20 
MeV to 300 GeV with an effective area peaking at ~ 9000 
cm 2 , giving it both detection area and energy maximum 
about 10 times those of EGRET. Fermi has been operated in 
survey mode, in which it views the entire sky every 3 hours, 
since shortly after its launch, and will continue to be in this 
mode for most of its > 5-year lifespan. The MAGIC experi- 
ment was recently upgraded to its second phase, MAGIC-II, 
and now consists of 2 large IACTs, each with a mirror area 
236 m 2 (the effective detection area for gamma rays is much 
larger, > 10 5 m 2 at optimal energies). The MAGIC tele- 
scopes are designed to have a low energy trigger threshold, 
<50 GeV near zenith, and are capabl e of repositioning to 
any point on the s ky within 30 seconds (jBastieri et alj|2005t 
I Albert et al1l2007T ). In this work, we will model the obser- 
vational capabilities of stereoscopic observations with this 
telescope using two different assumptions about instrument 
sensitivity at low energies. 

In Section O we present a model for high energy emis- 
sion from GRBs based on Swift and CGRO observations. 
We also briefly review recent predictions for the cosmolog- 



ical gamma-ray opacity due to pair-production interactions 
with UV and optical background photons, and the proper- 
ties of the instruments for which we are predicting detection 
rates. Our main results are presented in Section[3] where we 
show annual predictions for high-energy gamma-ray counts 
from GRBs for Fermi (Section and MAGIC (Section 
1 3 - 2 p . In Section[4] the impact of our findings for future obser- 
vations and EBL constraints is summarized. We also review 
our model in the context of the bright Fermi GRB 080916C, 
and calculate the number of counts that could have been 
seen by the MAGIC telescope had this burst fallen within 
its field of view. 



2 MODEL 

In order to estimate the number of GeV gamma-rays from 
GRBs that will be available to Fermi and MAGIC, we de- 
velop a simple model to estimate the fluence that could be 
seen by these experiments over a given time. As only bursts 
with known redshift are useful to our ultimate goal of prob- 
ing UV and optical background fields via attenuation effects, 
we base our analysis on the population of bursts observed by 
the Swift Burst Alert Telescope (BAT) with measured red- 
shift. Data for these bursts has been taken from the Swift 
GRB table available online Q. Our model rests upon two 
assumptions about high-energy emission: (1) that the pop- 
ulation of bursts seen by Fermi and MAGIC with redshifts 
that are eventually determined will be similar to these Swift 
bursts in number and fluence statistics, and (2) that these 
bursts produce high energy emission, both prompt and after- 
glow, with a fluence that is proportional to that observed in 
the BAT energy range, 15-150 keV. Here we define prompt 
emission as being approximately synchronous with the bulk 
of the flux seen by the BAT (T90, as discussed below), and 
the afterglow as occuring in the minutes to hours afterward. 
Our model is purely observational and phenomenological, 
and does not attempt to quantify in any way the intrinsic 
parameters of the bursts, nor do we make assumptions about 
the actual population statistics of these events. 



2.1 GRB Emission 

The popula tions of bu rsts seen by BATSE and Swift were 
analyzed in iDail l|2009h . who argued that these populations 
were similar, and that there was no difference between the 
subsets of optically detected Swift bursts with and without 
redshifts. This suggests that the events for which we now 
have redshift information are not different from the GRB 
population as a whole. We necessari ly exclude from our anal- 
ysis t he population of 'dark bursts' ( Virgili. Liang, fc Zhand 
l2009h . for which it is not possible to measure redshifts due 
to lack of an optical afterglow. We note that there have been 
suggestions that the low luminosity popu lation of bursts ar e 
distinct from their brighter counterparts (|Liang et al.ll2007h . 
The BATSE sample should be very similar to that viewed by 
GBM on Fermi, and this should enable a test of assumption 
(I)- 

We use the 145 bursts that were observed by Swift BAT 
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between January 2005 and June 2009 and have known red- 
shift. Figure [1] shows these bursts plotted as a function of 
redshift and BAT fluence. A considerable amount of fluence 
from these bursts arises largely from a few bright events; 
the brightest 10 per cent of bursts in the sample accounts 
for approximately 55 per cent of the fluence. While high en- 
ergy flux has only been seen from a handful of bright GRBs 
using EGRET, the fact that these bursts account for a large 
fraction of fluence seen at lower energies means that our 
assumption (2) should be reasonable even if the proportion- 
ality does not hold for faint bursts. We have not included 
Fermi LAT bursts such as GRB080916C in this analysis, 
although we do show where this event would have been in 
Fig. [1] based on its GBM fluence. We will discuss this event 
in the context of our emission model in Section 14.11 

Two recent papers have related the keV/MeV flux from 
GRBs to high energy emission, and have estimated t he ra- 
tio of fluences in these regimes. iLe fc Dermer! (|2009h esti- 
mated the count rate for GeV photons in the Fermi LAT 
based on the bursts seen by the EGRET spark cham- 
ber. The Flat / Fbatse fluence ratio inferred there var- 
ied from 0.05 to over 0.3. Based on the deadtime factors 
affecting some EGRET GRB observations these authors 
argue that a ratio of greater than 30 per cent between 
BATS E and EGRET is reasonable. lAndo. Nakar. fc Saril 
(2008) made the assumption that there is a log-normal dis- 
tribution of FGev/F Me v « Fegret / Fbatse in the roughly 
100 BATSE bursts that were in the field of view of EGRET. 
A maximum likelihood fit to the available data suggested 
a ratio of 0.003 ^ F G ev/F Me v «S 0.06. Slightly different 
assumptions about the high-energy spectrum and energy 
range of the GeV emission were used in each case. Le & 
Dermer used an index of -2 , similar to th e best fit to the 
EGRET bursts of -1.95 from lDingusl (|l995l ). while the Ando 
et al. work assumed an spectral index of -2.4 in the EGRET 
energy range. The latter noted that hardening the spectral 
index could have increased their values for the flux ratio 
coefficient significantly. 

In this paper we use the energy fluences seen by Swift 
BAT to predict GeV emission using the ratio 



p = -Fegret/Jbatse 



(1) 



where Fegret and Fbatse refer to the fluence (time- 
integrated flux over the duration of the event) across the 
EGRET spark chamber and BATSE energy ranges, taken 
to be 100 MeV - 5 GeV and 20 keV - 2 MeV respectively. 
A constant spectral index is assumed to be valid from the 
EGRET energy range up to > 100 GeV. We take an ap- 
proach similar to Le & Dermer and use a value of p = 0.1 
for prompt phase emission, which does a reasonable job 
matching the high energy fluence seen for the recent GRB 
080916C, see Section 14.11 This is higher than the range of 
values for p proposed by Ando et al., but that paper also 
assumed a softer high energy spectrum in deriving results. 

Afterglows are also a possible source of high-energy 
emission, though one that is even more poorly constrained 
than the prompt phase. There are various mechanisms that 
have been hypothesized as possible sources of GeV photons. 
A popular assumption invokes inverse-Compton upscatter- 
ing of synchrotron photons in the GRB outflow (SSC mech- 
anism) , although a variety of other sources are possible, such 
as SSC emission from the internal x-ray flares seen in after- 



glows or Compton upscattering of t hese photons by electrons 
accelerated in the external shock (|Fan et al.ll2008l ). Limits 
from EGRET observations suggest a typical flu ence ratio of 
0.01 to 0.1, and a spectral index of -1.5 to -2 (|Ando et al.l 
2008). For afterglows, we assume a ratio of p = 0.01 in this 
work. 

In converting the fluence seen by Swift BAT (15 - 
150 keV) to BAT SE flux, we assume a common Band 
jBand et all Il993l ) functional form over the BAT-BATSE 
energy range: 



dN 
dE 



A [E^ e" 
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(2) 



Here at\ and Q2 are the low and high energy indices, Ei, r is 
the break energy, and Q is the Heaviside step function. Un- 
fortunately, the relatively narrow energy band of the BAT 
does not allow one to resolve the structure of the Band peak 
in most cases for the Swift sample. We use -1 and -2.2 for the 
low and high indices, and assume a break energy of 250 keV; 
these are the typical values seen in an analysis of BATSE 
bursts by iPreece et all l|2000D . and lead to a ratio of 4.6 be- 
tween the BATSE and BAT fluences. The high energy flux 
in this model, effective in the EGRET energy range and 
at GeV energies, is taken to be a power law with a nor- 
malization found from the MeV-GeV proportionality, and 
is separate from the high-energy Band slope 02- The GeV 
flux is then 
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where f3 is the high-energy spectral index and is determined 
independently of the Band function parameters. We assume 
a high-energy spectral index of —1.95 for prompt phase pho- 
tons, consistent with the EGRET results, and a harder spec- 
tral index of —1.5 for the afterglow component, which would 
occur if the spectral peak from inverse Compton emission is 
at energies equal to or higher than those being observed. 
We do not discriminate between long and short bursts in 
our analysis. As can be seen in Figure [2] only a small frac- 
tion of the bursts in our sample can be classified as short 
(Tqq ^ 2 sec), and their contribution to the total flux is very 
low. 

The inverse-Compton scattering of photons to high en- 
ergies is limited by Klein-Nishina suppression, which reduces 
the flux of photons at energies that are higher than the elec- 
tron rest mass m e in the particle's rest frame. We will not 
include a possible cutoff due to this effect in our calcula- 
tion, but will reserve discussions of th e implications un- 
til Section 14.21 In the simple SSC case (Ch iang fc Dermerl 
Il999l : IZhang fc Meszarosll200ll : ISari fc Esinll200lh . this can 
be written as 



Ekn ~ Tb-y e m e c 



(4) 



where Ft is the bulk Lorentz factor of the outflow, and y e 
is the typical gamma factor for the electrons responsible for 
the synchrotron peak. Including the effect of redshift, the 
obse rved gamma-ray energies affected by this suppression 
are (Panaitcscu 200j|) 
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Figure 1. The bursts seen by Swift BAT for which we have well-measured redshifts, shown on axes of redshift and flucnce. The fluences 
observed by BAT at energies from 15—150 keV have been converted to the BATSE energy range (20 keV - 2 MeV) using a Band function 
over these energy ranges with peak 250 keV and low and high energy spectral indices of -1 and -2.2, respectively. For reference, we also 
show the position of GRB 080916C on this plot based on its G BM fluence (single red star). This burst was not seen by Swift until nearly 
17 hours after its detection by GBM and LAT ((Kennea 200^), and was not included in our analysis. 



Constraints from beaming and escape of high-energy radia- 
tion suggest large Lorentz factors for the bulk flow of ma- 
terial in t he prompt and early afterglow phases of bursts, 
F b ~ 100 (|Meszarosi r2006'). The electrons will typically have 
a power law distribution in energy determined by the cool- 
ing rate and therefore Klein-Nishina effects do not lead to 
an abrupt spectral cutoff, but if the typical Lorentz factor 
is sufficiently low then observations in the 10-1 00 GeV en- 
ergy r ange could be impacted. As mentioned in lAndo et alj 
l|200ct ). the electron Lorentz factor in external shocks in the 
afterglow is expected to be higher than that of the prompt 
emission, and therefore this may be a more likely source of 
detectable high energy photons. 



2.2 Cosmological gamma-ray opacity 

In this work we consider the effect of the UV-optical back- 
ground on high-energy GRB observations, using three dif- 
ferent models of the e volving background. T hese models are 
based on the work of I Gilmore et al] l|2009l ), in which semi- 
analytic models of galaxy formation and evolution are com- 
bined with estimates of quasar emissivity, ionizing escape 
fraction from star-forming galaxies, and processing of ioniz- 



ing photons by the inter-galactic medium to accurately com- 
pute both ionizing and non-ionizing fluxes. The low model 
in this work was based upon a WMAP3 cosmology that 
leads to delayed structure formation, due primarily to a low 
value for the dark matter power spectrum normalization, 
as = 0.761. The star-formation history in this model peaks 
near redshift 2 and is in agreement with the lowest values 
determined for high-redshift star formation. The rapid de- 
cline in star-formation rate density in this model leads to 
a UV background that is consistent with Lya ionization 
measurements only with a high escape fraction of ioniz- 
ing photons from galaxies, indicating that this is the min- 
imal level of UV background which can be consistent with 
high redshift observations. The fiducial model is based upon 
a semi-analytic model that uses a concordance cosmology 
(fi A = 0.7, tt m = 0.3, a 8 = 0.9). This model produces a UV 
luminosity density in good agreement with measurements 
at low and high redshift, and a star- formation rate peaking 
at z ~ 2.75. Finally, the fiducial high-peaked ('high') model 
features a star formation rate density and UV emissivity 
that is higher than that of the fiducial model at z > 3. The 
star-formation history in this model peaks at z = 5, and is 
consistent with some of the higher determinations of star- 
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formation rate density above redshift 3. Because this model 
provides maximum stellar emission at high redshift, and pro- 
duces a stellar mass density exceeding most observations, 
this should be considered a somewhat extreme model re- 
sulting in maximal background levels. This model converges 
to the fiducial model at redshifts lower than 3, and opacities 
for gamma-ray sources closer than this are not significantly 
different from the fiducial case. 



2.3 Instrument Properties 

For our purposes the performance of the LAT instrument 
on Fermi can be described with a relatively small number 
of parameters. We take the LAT to have an effective area of 
9000 cm 2 up to an energy of 300 GeV. While this may be 
an overestimate of the effective area at the highest energies, 
the number of photons seen at these energies is likely to be 
negligible. The integrated field of view for Fermi is found 
to be approximately ~ 20500 cm 2 sr; we therefore assume 
a field of view 20500/9000 « 2.28 sr. It is conservatively 
assumed in our analysis that Fermi will be in survey mode 
at all times, and that triggered rotations to view GRBs will 
not significantly raise the number of high-energy photons 
gathered. 

The observations of GRBs by IACTs such as MAGIC 
are highly sensitive to the capabilities of the instrument. 
While these telescopes have much larger effective collection 
areas than space-based instruments such as the LAT, other 
constraints such as the energy threshold, duty cycle, and 
time to respond to an alert must be taken into account in 
the analysis. The much larger effective area of these tele- 
scopes compared to Fermi is compensated by the relatively 
small probability that any single event will be observable. 
The effective area of IACTs is a strong function of energy 
in the sub-TeV regime, as the combined effects of the trig- 
ger efficiency and the reconstruction efficiency of low-energy 
showers in the analysis leads to a sharp decrease in effec- 
tive area near threshold. The potential of detecting GRBs 
is strongly affected by the low energy capabilities of the in- 
strument, due to the rapidly increasing opacity of the uni- 
verse to gamma rays above a couple hundred GeV for all 
but the closest bursts. For MAGIC, we have used the 'after 
cuts' form fo r the eff ective area as a function of energy from 
I Albert et all (|2008bl ). which is defined for angles near zenith. 
For observations away from zenith, the telescope loses sensi- 
tivity at low energies due to the increasing amount of atmo- 
sphere through which the particle shower is being observed. 
To estimate the change in the effective area function as a 
function of angle from zenith 9, we use the following simple 
approximation: 



A eff (E,0) = A eff (E-E', = 0), 
where 

E' = 50 GeV ■ ((cos 0)" 2 ' 5 - 1). 



(6) 



(7) 



This expression is based on a fit to the energy threshold 
re sults a t lowe r angles from zenith (8 £ 40 deg) presented 
in iFirpol l|2007l ). Additionally, we implement in our analysis 
an absolute cut at energy Emin, below which the sensitivity 
is zero and no gamma rays will be observed. We will per- 
form our analysis for values of 50 and 100 GeV in E m i n , to 



allow for uncertainty in low-energy sensitivity and data re- 
construction. Note that this parameter is distinct from the 
instrument energy threshold in that it is an absolute cut and 
is not dependent on the spectral properties of the source. Re- 
cent esti mates of the instr ument sensitivity for stereo obser- 
vations (| Colin e t al. 2009) suggest that our more optimistic 
cut of 50 GeV may be reasonable for future observations. 

A realistic estimate of the instru ment duty cycle ident- 
ical for our analysis. As outlined in iBastieri et alj l|2005h . 
there are several requirements for the operation of MAGIC, 
including distance of the sun from zenith (> 108 deg), a 
minimum angular distance of the moon from the observa- 
tion field (> 30 deg), and humidity and wind requirements. 
For the duty cycle of the instrument, we use the standard 
value of 10 per cent, and note that this is supported by the 
fraction of GRBs (9.2 per cent) which we re observed during 
follow -up in 42 months of observations (|Garczarczvk et aU 
120081) , and is similar to reports from observations with other 
IACTs. 

A major challenge for ground-based attempts to detect 
GRBs has been the highly transitory nature of the emis- 
sion. Responding to an event requires the minimization of 
the individual components that contribute to the total delay 
time, including the time for the detecting satellite to confirm 
a burst in progress, the time to transmit this information, 
and then the time for the ground-based telescope to slew 
to the coordinates and begin taking data. The first quan- 
tity depends on the satellite response time and strength of 
the burst. Typical numbers for the Swift BAT are < 15 sec 
to confirm and transmit a coordinate with a precision of a 
few arcminutes. Communications of these c oordinates are 
receiv ed by the IACTs in real time (~ 2 sec, IBastieri et al.l 
120051) over the GRB Coordinate Network (GCNfl However, 
the final step of repositioning a telescope the size of MAGIC 
on the time scale required to view prompt GRB emission 
presents a major engineering challenge. Ensuring personnel 
safety is also a major practical concern in minimizing re- 
sponse time, which requires that the telescope be able to 
reposition without warning at any time during operation. 

For our analysis, we assume a delay time Tdeiay in 
prompt observations, which incorporates all three times dis- 
cussed above. We assume a typical report time of 15 seconds 
for the burst alert to reach the instrument, and, in the case 
of MAGIC, a 30 sec slew time to move to the target and be- 
gin observations. The total of 45 seconds is about equal to 
the lowest time, 43 sec, reported in lGarczarczvk et al.l (|2008T l 
for all of MAGIC GRB responses to date, and is therefore 
optimistic. To compute the flux in the prompt phase of a 
GRB that can be seen by ground-based observations after 
this delay time we use the T90 variable in BAT flux reported 
for our sample of bursts, which is the time elapsed between 
the arrival of the first 5 and first 95 percent of the total 
GRB fluence in the BAT band. The fluence of the burst is 
modified by a factor 



F = Fbat ■ MAX 



Tgo — T c 



delay 



(8) 



that is, we take the prompt phase emission profile to be ap- 
proximately flat over this timescale, and reduce the fluence 
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Figure 2. The fiuence factors arising due to the slew time delay 
for ground-based telescopes. The solid black curve shows the per- 
centage of prompt fiuence in the total sample which is seen after 
a time T^ e i ay , if the GeV flux is proportional to that in the BAT 
energy range and is constant across Tgo. The red curve shows the 
fraction of GRBs for which Tgo is longer than Td e i a y This plot 
shows the percentage of high energy emission that will be missed 
by a telescope with a given delay time (due to slewing and other 
factors) in our model, after averaging over all Swift bursts. 

by the proportion of the prompt phase that was missed. De- 
lays due to IACT response time do not affect the observation 
of afterglows. We also do not account for the fact that on a 
timescale of hours the rotation of the Earth can bring after- 
glows into and out of the viewing region of the telescope. In 
Figure [2] we show how the time delay in IACT observations 
affects observed flux in our model. It is worth pointing out 
that we do not find Tgo to be correlated with either fiuence 
or redshift in our sample; its correlation coefficients with 
these variables are -0.017 and -0.014, respectively. 



3 RESULTS 

Based on the model for GRB emission and instrument per- 
formance developed in the previous section, we predict the 
number of high-energy gamma rays from GRBs that will be 
seen by Fermi and MAGIC, using the data set of Swift BAT 
GRBs with confirmed redshifts and the emission model dis- 
cussed in the previous section. In Table [T] we review the 
parameters for the emission model and instrument proper- 
ties that we are using in this section. 

3.1 Predicted Fluences for Fermi 

With its wide field of view, Fermi is expected to view several 
high-redshift GRBs per year. We therefore present results for 
this section in terms of broad bins in redshift, so as to have 
a reasonable statistical sample in each bin. In Table [2] we 



Table 1. Some of the parameters we use for calculations of this 
section. See Sections 12. II and 12. 31 for more details. 



Fbatse/-Fbat 4.6 BAT to BATSE fiuence ratio 

p pr 0.1 Fiuence ratio for prompt phase 

f3 pr -1.95 Prompt phase spectral index 

p ag 0.01 Fiuence ratio for afterglow phase 

(Sag -1.5 Afterglow spectral index 

Tdeiay 45 sec Delay time for MAGIC observations 



Table 2 . The redshift bins we use in our analysis and the numbers 
of GRBs and total fiuence in each bin for the sample of Swift 
bursts. The fluences shown are in the BAT energy range, 15-150 
keV. These data are over 54 months of observations, from January 
2005 to June 2009. 



Redshift Bin 


N(GRB) 


Fiuence (10 7 erg cm 2 ) 


1 < z < 2 


29 


1595.9 


2 < z < 3 


33 


864.9 


3 < 2 < 4 


21 


665.1 


4 < z < 6 


10 


228.3 



show the redshift bins, number of bursts seen in 54 months 
operations, and total fiuence. 

After calculating the high energy gamma-ray fiuence 
for each burst in the sample, the flux is attenuated us- 
ing optical depths calculated fro m the evolving back g round 
sp ectral energy d i stribu tion; see iMadau fc Phinnevl (|l996l ) 
or iGilmore et all (|2009h for a review of these calculations. 
In each redshift bin, we use attenuation factors averaged in 
redshift, 

~dE~ = Tie J zi T^T dz > (9) 

where dF/dE is the high energy spectrum from Equation [3] 
and t(E, z) is the EBL optical depth as a function of ob- 
served gamma-ray energy and source redshift. High-energy 
gamma rays are therefore assumed to originate from sources 
evenly distributed across the redshift bin. 

Figure [3] presents the main results of this section; the 
mean number of photons above a given energy available 
without absorption by the EBL, and after absorption by our 
background models. These predictions are made by combin- 
ing our models for high-energy GRB emission for the Swift 
population and Fermi instrument properties, making adjust- 
ments for the field of view and effective areas for the two de- 
tectors. Our results are divided into the four redshift bins, 
and show the averaged total amount of fiuence from sources 
at these redshifts per year. 

In Figure U we show how the mean fiuence predictions 
can be expected to vary from year to year, based on the 
number of bursts in the Swift sample and the variance in 
fiuence in this population. This figure has been created using 
a simulation of LAT observations, and assuming randomly 
occurring bursts with the flux and redshift distribution of 
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Figure 3. The mean number of integrated photons above a given energy visible to Fermi per year, up to the maximum LAT energy of 
300 GeV. The dotted line shows the unabsorbed rate, whi le the solid blue and black (upper and lower, respectively) lines show results 
with attenuation due to the low and fiducial EBL models of lGilmore et al.l ll2009f) . The solid orange line (lowest line in the bottom panels) 
shows the integrated counts using the high model, which leads to greater opacity than the fiducial model above redshift 3. 



the Swift population. The upper two distributions show the 
number of high-redshift (1 < z < 6) bursts falling in the field 
of view of the detector over a period of 1 (upper-left) and 5 
(upper-right) years; this follows a Poisson distribution. As 
before, we do not account for the possibility of the spacecraft 
autonomously slewing to view events with the LAT after 
being triggered by the GBM or another experiment. The 
lower panels show how the stacked fluence collected can be 
expected to vary from the predictions in the previous plots. 



3.2 Predicted Fluences for MAGIC 

For MAGIC, we begin by considering the year-to-year prob- 
ability that a given number of high-redshift GRBs will oc- 
cur in a region of sky where they can observed with a low 
energy threshold. By multiplying the duty cycle of the in- 
strument (10 per cent) with the sky coverage (11.7 per cent 
for 8 m ax = 40 deg), we find that only ~1 per cent of bursts 
can be observed with a reasonably low energy threshold. 



Our sample consists of 145 bursts seen over 54 months, 96 
of which are at z > 1, and therefore the expected number 
of bursts per year is somewhat less than 1. Therefore, we 
begin by predicting the probability that any bursts will be 
visible. Once we understand this probability, we will look at 
the photon statistics for a single burst. 

The sky coverage of the telescope increases approxi- 
mately as the square of the maximum allowed angle from 
zenith. However, we find that the number of photons pre- 
dicted from distant [z > 1) bursts does not increase signif- 
icantly beyond an angle of about 40 degrees, as the energy 
threshold of the instrument rises above the energies at which 
the universe is transparent to gamma-rays. GRBs which oc- 
cur far from zenith are therefore shrouded from view by EBL 
attenuation unless they are at low redshift. The number of 
bursts per year will also depend on whether prompt or af- 
terglow phases are being considered, as the Tgo duration of 
some GRBs will be less than Td e i ay , preventing most or all of 
their prompt emission being viewed. In Table [3] we present 
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Figure 4. The variation in yearly predictions for observations with the Fermi LAT of GeV gamma rays from GRBs with determined 
redshifts. The upper-left plot shows the annual probability of a given number of GRBs occurring in view of the LAT (taken here to be 
angles $J 50 deg from boresight) with redshifts between 1 and 6. This does not take into account the expected number of GeV photons 
from each GRB, which may be less than one. The lower-left plot is the yearly distribution of stacked fiuences, normalized to the mean 
yearly predictions in Figure [3] On the right-hand side we show the same quantities computed for a 5— year period. 



Table 3. The annual probabilities of a given number of GRBs falling within a given angle from zenith, assuming a duty cycle of 10 
per cent. Under 'afterglows', we include all bursts, and for 'prompt' only those with Tgo > Tdelay Note that this calculation is for the 
number of GRBs falling in view of the instrument, and does not predict the number of expected photons, which may be less than one. 
In Figs. [6] through [9] we show the distribution function for the number of photons received from a single burst falling in view of the 
instrument, using our emission model and taking into account EBL attenuation and instrument effective area. These factors also take 
into account the antisolar bias found in the distribution of Swift GRBs. 







Prompt {T de i ay = 45s) 






Afterglow (T delay = 0) 




N(GRB) yr" 1 


20 deg 


30 deg 


40 deg 


20 deg 


30 deg 


40 deg 





0.95 


0.90 


0.82 


0.92 


0.82 


0.71 


1 


0.048 


0.099 


0.16 


0.079 


0.16 


0.24 


^2 


0.0012 


0.0057 


0.016 


0.0036 


0.016 


0.046 
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Figure 5. The distribution of the angular separation 9 S e P be- 
tween Swift GRBs with redshifts and the position of the sun. 
This is plotted as cos(8 se p), which would yield a flat distribution 
if the GRBs appeared independently of solar position. The dot- 
dashed line is the fit used in determining the enhancement factor 
for IACT observations due to the solar anti-bias in burst alerts; 
see text for details. 



0.30 - 




counts 

Figure 6. Predictions for the number of gamma-rays visible to 
MAGIC from a high-redshift (z > 1) GRB; see Table [3] for the 
rates at which these events are expected to occur in view of the 
telescope. Here we use parameters T £ j e ; aH =45 sec and 9 max = 40 
deg, and assume minimum energy of 50 GeV in telescope sensitiv- 
ity. Blue, black and orange lines show the results after attenuation 
by the low, fiducial, and high EBL models, respectively. The left- 
most bin is the probability of zero photons being received. 

probabilities for the number of GRBs visible to MAGIC per 
year, within a given maximum angle from zenith. 

Another factor that must be taken into account for this 
analysis is the anti-solar bias in the distribution of Swift 
GRBs. As the determination of redshifts from afterglow ob- 
servation is hampered by glare from the sun, Swift pref- 
erentially finds GRB events in the anti-solar direction. This 
works to the advantage of IACTs, which can only operate at 
night, by effectively increasing the duty cycle with respect 
to Swift-triggered bursts. As low-threshold IACT observa- 
tions are limited almost exclusively to solar angles of >90°, 
we can estimate the effect of this bias factor on the duty 
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counts 

Figure 7. As in the previous figure, but with E m i n = 100 GeV 

cycle by taking the GRB distribution to be flat for separa- 
tion angles >90° and modeling a cutoff in GRB detection 
rates at smaller angular separations. In Fig. [5] we show the 
distribution of bursts on the sky with respect to the position 
of the sun. Using the simple fit shown, we find an enhance- 
ment factor of ~ 1.38 for IACTs relative to an isotropic 
distribution. 

We find that the probability of a high redshift burst 
falling within the field of view of an IACT in any particu- 
lar year is small, ~30 per cent for 6 max — 40 deg. At this 
angle, the energy threshold is approximatel y double th at at 
zenith, and rises rapidly for larger ang les l|Firpoll2007l ). It 
is not realistic to expect to see GRBs at z > 1 at larger 
zenith angles for either of our EBL models due to attenua- 
tion, although low redshift bursts could be visible if photon 
emission takes place at high enough energies. In years in 
which a nonzero number of bursts is in view of the tele- 
scope, the predicted flux is expected to vary highly due to 
the large range of fluences seen in the Swift sample (y-axis 
in Figure [TJ. It is therefore not particularly useful to de- 
scribe an 'average' year, as we did in the previous section 
when discussing Fermi, as the median number of observed 
gamma-rays from GRBs in this case is zero. 

In Figures|6]and[7]we show as a histogram the probabil- 
ity distribution for the number of detected photons expected 
from a single distant (z ^ 1) GRB within 9 max - This calcu- 
lation has been done by repeatedly sampling our catalogue 
of Swift bursts, randomizing the sky position of the GRB 
within the disk of radius 8 m ax each time. The MAGIC ef- 
fective area as a function of energy is then calculated using 
Equation [6] We apply absolute cuts in energy of E c = 50 
and E c = 100 GeV to our analysis, respectively in these 
two figures. The predicted high energy flux is attenuated 
using our three evolving background models. Additionally, 
the redshift of each burst is randomized each time it is sam- 
pled by an amount —0.2 Az ^ 0.2. Figures [8] and [9] shows 
the same plots including only the prompt phase photons. 

Figure [TU] shows the expected dN/dE spectrum of pho- 
tons arriving from GRBs at different redshifts. This plot 
takes into account the assumed gamma-ray spectrum, the 
instrument effective area, and EBL attenuation for four dif- 
ferent redshifts. For MAGIC, we find that photons are ex- 
pected to be seen in a fairly narrow energy range peaking 
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Figure 8. Result for observations with an energy cut at E m i n = 
50 GeV, and here including only photons from the prompt phase 
of emission. Line types are as in Figure [6] 
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Figure 9. Result for observations with E m i n = 100 GeV, and 
including only photons from the prompt phase of emission. Line 
types are as in Figure [6] 

at ~ 100 GeV. The number of photons expected near the 
minimum energy of the instrument are suppressed due to the 
small effective area at these energies. At higher energies, the 
spectrum declines rapidly due to EBL attenuation. As none 
of these spectral factors depend on absolute fluence in our 
model, these results are valid regardless of the luminosity of 
the GRB. We have normalized the spectra to the average 
yearly flux for convenience, however this normalization in 
itself has little meaning due to the large amount of variance 
in our predictions. 



4 DISCUSSION 

We have attempted to make predictions for the number of 
high-energy gamma rays that can be seen by current gen- 
eration of telescopes targeting the GeV energy range. Al- 
though we have strived to design a simple and straightfor- 
ward model, there are a large number of uncertainties in 
predicting the high-energy emission of GRBs. Over the next 
few years, observations in the GeV energy band will be able 



to constrain many of the assumptions that we have used 
here. 

Our findings suggest that the Fermi LAT will typically 
observe at least 3 to 4 bursts per year with redshifts that 
are determined to be greater than 1. Over the lifetime of the 
mission (>5 yr), this means that multiple events will likely 
be seen in each of the redshift bins. While not all these GRBs 
will necessarily have detectable GeV emission, the stacked 
results could yield enough photons above 10 GeV to place 
constraints on the EBL, and possibly differentiate between 
the two models we have presented here. It is difficult to make 
generalizations about yearly predictions, due to the large 
amount of variance in both number of bursts viewed and 
fluence per individual burst, which vary over orders of mag- 
nitude. Our predictions become more stable over the 5-year 
instrument lifetime (right-hand panels of Fig. [4| ; the proba- 
bility of our results varying by more than a factor of ~ 2 from 
predictions on this timescale is small, at least when totals 
from all redshift bins are considered (this is likely smaller 
than the expected variation in p from Section We pre- 
dict, on average, several photons per year above 10 GeV for 
1 ^ z ^ 2 and two or three from 2 ^ z ^ 3. This is com- 
parab le to the more optimistic predictions of iLe fc Dermerl 
(2009), who considered all long-duration bursts over the full 
sky with the same high energy spectrum and MeV-GeV flu- 
ence ratios. We have not included autonomous repoints of 
the instrument in response to a GBM trigger in our model. 
If repoints are performed frequently and prove an effective 
way to view GeV emission from bursts, then this could effec- 
tively increase the LAT field of view to 1/2 the sky, a factor 
of 2.5 higher than we consider. 

For an IACT like MAGIC, the annual likelihood of view- 
ing one or more distant GRBs within 40 degrees of zenith is 
less than 30 per cent. We have chosen this cutoff due to the 
rapidly growing energy threshold at higher angles, and the 
fact that the EBL absorption attenuates most emission at 
these energies for sources at z > 1 in all of our EBL models. 
As shown in Figure [101 the majority of photons from these 
bursts can be expected to arrive at ~ 100 GeV, or lower 
for high redshift sources. This means that it is the capa- 
bilities of the detector near to the lower end of its energy 
range that are most important to viewing distant GRBs, 
not the energies at which the instrument necessarily has the 
most sensitivity. As can be seen in Figs. [7] and [9] increas- 
ing the assumed minimum energy of the instrument to 100 
GeV greatly decreases sensitivity to GRBs, especially when 
combined with a higher realization of the EBL. Sensitivity 
below 100 GeV is therefore critical to detecting high redshift 
GRBs if the UV-optical background is as high or higher than 
found in our fiducial model. 

In this analysis we have restricted ourselves to GRBs 
at z ^ 1. For GRBs at lower redshifts, more photons may 
be observable in the 100 GeV to 1 TeV decade, but these 
bursts only represent a minority of those for which we have 
redshifts. Another reason for focusing on the 10 to ~100 
GeV decade is that we are relying on EGRET observations 
at maximum energies of ~ 10 GeV to model our high-energy 
emission, and we expect this model to become increasingly 
uncertain at higher energies. Our main results from this sec- 
tion show the potential of receiving a large number of photon 
counts from a single GRB. Even a relatively small number 
of gamma rays could be very useful in constraining the EBL 
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Figure 10. Predictions for the spectrum observable by MAGIC from distant GRBs, at 4 different rcdshifts. Here we assume observations 
at 20 degrees from zenith with an E m i n = 50 GeV. The dotted black line is the unattenuated spectrum after including instrument 
effective area, and the solid orange, black, and blue lines (lower to upper) show the flux after attenuation by the high, fiducial, and 
low EBL models. The fluence is normalized to the average per year, as was done for our Fermi results. However, as we have seen, the 
year-to-year prediction is highly variable, so the vertical scale of the plot should be taken as arbitrary, and Figs. [6] through [9] used as a 
gauge of expected number of counts. 



through its effects on the spectrum and total power of the 
VHE emission. Our results suggest that there is a large de- 
gree of chance involved in seeing GRBs from the ground, but 
there could be a large payoff for our knowledge of cosmology 
from even one success. 

For both MAGIC and Fermi, we have based our calcu- 
lation on the population of bursts seen by Swift for which 
redshifts were eventually determined. However, Fermi is a 
capable finder of transients, covering approximately 1/5 
of the sky with the LAT and possibly more after GBM- 
triggered repointings are taken into account. Above 10 GeV, 
the LAT has an angular resolution of sj 0.1 degrees, allow- 
ing strong bursts to be targeted for multiwavelength obser- 
vations. It is therefore possible that Fermi will enable the 
determination of many GRB redshifts without simultaneous 
Swift detections, and the number of high-redshift bursts we 
have predicted for the next few years could be an underes- 



timate. This could improve the prospects for detecting GeV 
photons from these sources with both Fermi and IACTs. 



4.1 Simulated Results for GRB 080916C 

One exciting implication in our findings is the potential pay- 
off from a single bright GRB. For MAGIC, we find that while 
the probability of seeing photons from any single event is 
quite small, the reward for catching a burst of intermediate 
to high fluence could be hundreds or thousands of photons 
observed within a narrow energy range. The detailed spec- 
trum from such an event could be invaluable for constrain- 
ing the UV background and high-redshift galaxy formation. 
Just as a demonstration of how many photons one event 
could provide, we consider recent GRB 080916C which was 
observed in its full prompt phase by the Fermi LAT, as 
well as GBM. Swift did not observe this event until nearly 
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a day afterwards, and it was not included in the calcula- 
tions of the previous section. This burst, which occurred on 
September 16, 2008, is among the brightest GRBs ever seen, 
and with a redshift of 4.35 ± 0.15 it is the most energetic 
b urst currently know n l|Greiner et al.| [2009). As described 
in lAbdo et al.l l|2009l ) , it was seen by the L AT at an angle 
of 48 degrees from boresight following a trigger from the 
GBM. A total of 145 gamma rays above 100 MeV and 14 
gamma rays above 1 GeV were reported. The last of those 
14 gamma rays arrived approximately 46 seconds after the 
initial trigger. The highest energy gamma ray was measured 
to be 13.22^j'g4 GeV, and occurred 16.54 seconds after the 
trigger. 

A useful test of our emission model is to ask if we ar- 
rive at similar results for the number of LAT-observed pho- 
tons above 1 GeV, using the GRB 080916C GBM fluence of 
1.1 x 10~ 4 erg cm -2 as the basis for our calculation. Follow- 
ing the same analysis as in Section f3. II for this single GRB, 
we predict 24 and 23 photons above 1 GeV respectively for 
our low and high EBL models, compared to 14 observed by 
the LAT. For photons above 10 GeV, we predict 1.9 and 
1.4 photons for the two models, compared to one seen. The 
number of predicted cumulative counts falls below 1 for en- 
ergies of 15 and 11 GeV. Thus, while our emission model 
overpredicts the number of photons around 1 GeV, it does 
a reasonably good job of predicting the highest energy pho- 
ton seen by the LAT for this event, though it is difficult to 
draw conclusions with such small statistics. The role of at- 
tenuation by the EBL in this energy range is minimal. Our 
low and high models give optical depths of 0.06 and 0.22 for 
a photon of 13.2 GeV from redshift 4.35, corresponding to 
attenuation factors (e~ T ) 0.94 and 0.80. 

Next, we consider a hypothetical observation of this 
GRB by MAGIC. As the last photon from GRB 080916C 
above 1 GeV arrived at about 46 seconds after the initial 
trigger, it is unlikely that MAGIC or other IACTs could 
have seen much of the prompt emission, and the counts we 
calculate here are entirely from a hypothesized afterglow 
component. For an IACT, observations of this burst would 
have been heavily dependent on the assumed EBL model, 
due to its high redshift. We calculate that with a 50 GeV 
energy cut at zenith MAGIC would have seen 350 gamma 
rays for the EBL in our low model; in the fiducial and high 
models 58 and 19 gamma rays would be seen. This assumes 
ideal viewing conditions, with the GRB occurring directly 
overhead, and in reality the chances of such an occurrence 
are exceedingly small. At higher zenith angles, the num- 
ber of observable photons declines rapidly due to the higher 
threshold and the EBL attenuation being a strong function 
of energy. If the event were instead seen at 45 degrees from 
zenith the predicted gamma-rays counts for those same three 
EBL models would be 61, 2, and 0.24. As we see, predictions 
can vary enormously for high-redshift GRBs depending on 
the background model. The fiducial and high models create 
a dense background of UV photons due to earlier star for- 
mation, and 100 GeV gamma-rays from a source at z ~ 4 
are attenuated by a factor > 100. Our low model has much 
less star formation at high redshift, and the optical depth 
to gamma rays is much lower, although the universe is still 
optically thick (r > 1) at this redshift for photons above 50 
GeV. 

As mentioned in lAbdo et al.l l|2009h . the high energy 



emission measured by the LAT was delayed slightly com- 
pared to the GBM flux. The highest energy photon, and 
two others that had energies above 6 GeV, did not arrive 
until more than 83 per cent of the prompt GMB fluence had 
been received (as seen in Table 1 of this reference) . While it 
is difficult to draw conclusions from one event, this may in- 
dicate that the VHE photons produced in the prompt phase 
may arrive later than the lower energy fluence which defines 
T90, or possibly that the spectrum hardens with time and 
GeV photons tend to arrive later than lower energy emis- 
sion. We have assumed in our analysis that MeV and GeV 
prompt-phase flux are directly proportional in time with a 
constant spectrum. If there is a delay or spectral hardening 
it could work to the advantage of ground-based instruments, 
allowing them more time to react to a GRB report than we 
have assumed here. 



4.2 Intrinsic Spectral Cutoffs 

One hurdle in detecting gamma-ray attenuation features 
could be the existence of a spectral cutoff due to either the 
Klein-Nishina cutoff or internal absorption of gamma rays. 
As described in Equation [5] the relevant energy scale is de- 
termined by the bulk Lorentz factor of the GRB ejecta and 
the typical electron relativistic Lorentz factor. The analysis 
of GRB 080916C suggests a bulk factor of T bu i k ^ 887 ± 21 
during the time intervals when the highest energy gamma 
rays were emitted. If the electron Lorentz factor was at least 
~ 10 3 , then emission that could be observed by MAGIC 
would not be affected. However, as the most powerful GRB 
on record, parameters for GRB 0809 16C may not be repre- 
sentative of the total sample. One potential danger is that 
the typical energy of the cutoff could exist at roughly the 
same GeV energies where we expect EBL attenuation fea- 
tures to be seen. Not only could a sharp spectral cutoff be 
mistaken for attenuation by background radiation, but the 
factor of (1 + z)" 1 from cosmology could mimic the redshift 
evolution of EBL attenuation. 



4.3 Future Experiments 

One reason we have restricted ourselves to current exper- 
iments in this discussion is that, as we have seen, the de- 
tails of instrument capabilities can have a large impact on 
predictions, and our results are most meaningful when we 
can incorporate well-tested and verified instrument param- 
eters into our model. But as the GeV emission of GRBs 
and subsequent extinction are certainly not questions that 
are going to be decisively answered by the current genera- 
tion of instruments, our discussion would not be complete 
without mentioning briefly a few important upcoming ex- 
periments. The next phase of the H.E.S.S. array will fea- 
ture a 600 m 2 mirror at the center of its current 4-telescope 
configuration; this central 'T5' telescope will be the largest 
IACT yet built. This upgrade is scheduled for completion 
later this year, and will lower the en ergy threshold down t o 
~30 GeV at zenith angle 18 degrees (|Becherini et al.l l2008). 
Over the next decade, several ground-based experiments 
can provide more sen sitivity to VHE photons from GRBs 
(Willi ams et al.ll2009Tl. The Advanced Gamma- Ray Imaging 
System (AGIS) ilBucklev et all 120081) and Cherenkov Tele- 
scope Array (CTA) ( Martinez 20081 ) are two future concepts 
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for IACT arrays that may be constructed during the next 
decade. Both of these arrays, when fully constructed, would 
have much larger collection areas than any current experi- 
ment, and would likely have energy coverage over most of 
the 10 to 100 GeV decade. The low energy threshold, < 20 
GeV, and greater sensitivity will enable detection of GRBs 
and measurement of attenuation by the EBL out to very 
high redshift, which could clarify the nature of high-redshift 
star formation. Unfortunately, these telescopes will not be 
able to overcome the intrinsic difficulties of the Cherenkov 
technique, namely low duty cycle, loss of sensitivity away 
from zenith, and the need to be triggered for transient ob- 
servations by another experiment. Our results suggest that 
due to the stochastic nature of GRBs, persistence may ulti- 
mately be the key to detecting these events from the ground. 
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